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I.  SUMMARY 

Magnetoplasmadynamic  (MPD)  arc jet  thrusters  inherently 
involve  close  coupling  of  the  electromagnetic  discharge  structure 
and  the  flow  field  in  the  thrust  chamber.  The  discharge  and 
internal  flow  field  structures  are  not  simply  determined  by  the 
arc jet  geometry,  but  depend  on  the  values  and  partitioning  of 
injected  mass  flow  for  each  operating  current  condition.  In 
order  to  understand  the  internal  flow  dynamics  of  an  MPD  arcjet, 
experiments  are  performed  in  which  electromagnetic  fields  and 
flow  properties  are  measured  within  the  thrust  chamber.  Electric 
and  magnetic  field  probes  are  used  to  obtain  the  electromagnetic 
structure,  while  spectroscopic  techniques  are  used  to  estimate 
particle  velocities  and  densities. 

Results  from  the  past  year’s  effort  include  measurements  of 
Doppler  widths  in  argon  and  hydrogen  as  functions  of  position  in 
the  thruster  for  two  arrangements  of  mass  injection:  100  %  at 
backplate  mid-radius  and  50-50%  at  backplate  mid-radius  and 
cathode  base.  Reduced  data  suggest  that  a  substantial  radial 
component  of  flow  exists  in  the  former  case,  while  little  radial 
flow  is  indicated  for  the  latter  case.  Radial  flow  could  result 
in  significantly  inceasing  shock  and/or  viscous  losses  near  the 
cathode.  Performance  improvements  in  MPD  arc jets  would  thus 
depend  on  properly  matching  the  thrust  chamber  injector  geometry 
to  the  flow  field  actually  obtained  with  particular  choices  of 
total  current  and  mass  flow. 
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II .  INTRODUCTION 

The  magnetoplasmadynamic  (MPD)  arc jet  has  attracted 
attention  for  many  years  because  of  its  promise  for  providing  a 
high  thrust  density,  high  specific  impulse  space  propulsion 
technique.  A  particularly  important  characteristic  of  the  MPD 
arcjet,  (in  contrast  to  electostatic  thrusters) ,  is  the 
simplicity  of  the  electrode/channel  structure  accelerating  and 
guiding  the  thruster  flow.  While  the  solid  structures  of  the  MPD 
thruster  can  be  quite  simple  (e.g.,  two  coaxial  electrodes),  the 
internal  flow  of  plasma  in  the  thrust  chamber  may  be  very  complex 
since  the  MPD  flow  inherently  depends  on  close  interaction  of  the 
electromagnetic  discharge  and  flow  structures.  These  flow 
structures  are  not  simply  determined  by  the  thrust  chamber 
geometry  but  can  vary  with  changes  in  overall  operating 
parameters  such  as  current  and  rate  of  mass  flow  injection.  In 
order  to  understand  the  performance  characteristics  of  MPD 
arcjets,  it  is  necessary  to  determine  the  electromagnetic 
discharge  and  flow  structures  within  the  thrust  chamber.  Such 
determination  requires  diagnosis  of  electromagnetic  field 
distributions  and  measurement  of  the  particle  densities  and 
velocities  in  the  thrust  chamber  flow. 

The  present  research  program  is  directed  toward  diagnosing 
the  MPD  thrust  chamber  flow  and  comparing  the  measured 
electromagnetic  and  flow  structures  with  theoretical  predictions 
(e.g.  MHD  code  results) .  The  electromagnetic  structure  is 
obtained  by  internal  electric  and  magnetic  probe  techniques.  Non- 
intrusive  optical  techniques,  particularly  time-and  spatially- 
resolved  spectroscopy,  provide  measures  of  the  particle 
velocities  and  chemical  states.  In  addition  to  experimental 
measurements,  a  set  of  calculations  using  the  MACH2  two- 
dimensional  MHD  code  [1]  has  been  started  based  on  the  actual 
arcjet  geometry  and  operating  parameters  of  the  experiment. 
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III.  PROGRESS 

In  the  past  year,  a  new  MPD  arc jet  system  was  Installed  and 
tested.  The  device,  shown  in  Figure  1,  has  a  brass  outer  anode 
8.41  cm  ID  and  a  copper-tungsten  alloy  cathode  3.81  cm  OD.  The 
channel  length  is  5  cm  from  the  boron  nitride  insulator  backplate 
to  the  exit  plane.  Argon  or  argon-hydrogen  mixtures  were 
delivered  to  the  arc jet  through  16  choked  orifices  in  the  boron 
nitride  backplate  at  the  mid-radius  of  the  channel  gap.  Current 
is  delivered  to  the  arc jet  from  the  PFN  with  a  pulse  width  of 
0.8  ms  and  peak  currents  of  up  to  45  kA. 

A  consistent  set  of  operating  conditions  were  used  for 
magnetic  probe  and  spectroscopic  measurements.  The  PFN  was 
charged  to  8  kV  and  the  gas  flow  rate  was  6  g/s  (argon  with  1.5 
percent  added  hydrogen) .  Figure  2  shows  the  current  waveform 
for  these  conditions.  The  quasi-steady  current  is  approximately 
20  kA.  The  arcjet  breech  voltage,  Figure  3,  attains  a  quasi¬ 
equilibrium  value  of  50  V  for  approximately  400  n s  about  200  ns 
after  current  initiation.  Two  magnetic  field  probes  were  used  to 
map  the  enclosed  current  contours.  Over  70  data  shots,  at  two 
data  points  per  shot,  were  taken  to  map  60  locations  in  the 
arcjet  thrust  chamber.  Figure  4  shows  the  enclosed  current 
contours  in  the  arcjet  operating  at  20  kA  and  6  g/s.  Note  that 
the  current  distribution  is  basically  retained  within  the  thrust 
chamber,  with  less  than  5  percent  of  the  total  current  beyond  the 
exit  plane.  There  are  also  indications  that  near  the  boron 
nitride  backplate  at  the  midchannel  gas  inlet  ports  the  current 
distribution  is  perturbed,  possibly  in  response  to  the  entry  of 
relatively  high  density,  low  speed  flow  into  the  thrust  chamber. 
The  J  x  B  forces  on  the  plasma  have  significant  radial  components 
near  the  anode  and  cathode  suggesting  radial  velocity  components 
and  indicating  possible  increases  in  plasma  density  towards  the 
conductor  surfaces,  particularly  the  cathode. 
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Figure  1.  Schematic  of  new  MPD  arc  jet  for  study  of  MPD  thrust 
chamber  flow  dynamics.  Gas  injection  is  through  16 
choked  orifices  at  mid-radius  in  boron  nitride 
backplate. 


4 


VOLTS  KILOAMPS 


20 


I 

I 


Figure 


TIME  (MILL ISECONOS) 

.  Current  delivered  to  the  MPD  arc jet  has  a  FWHM  of 

900  ms  with  a  risetime  of  100  m s.  An  initial  voltage 
of  8  kV  on  the  PFN  gives  a  current  of  approximately 
20  kA. 


80 

60 

40 

20 

O 

-20 

.  Breech  voltage  attains  a  quasi-equilibrium  value  for 
approximately  400  ms  about  200  ms  after  current 
initiation. 


0.0  0.2  0.4  0.6  0.8  1.0 

TIME  (MILLISECONDS) 


5 


I 


UJ 

a 

o 

X 


o 


UJ 

zq< 

Ohj 

ttcci 

9  id 

CDmO 

z  < 
m 


H  H 
UJ  2 
WJH 

<zo 
ts  m  a. 


© 

jC  iH 
4->  >0 
-H 
P  TJ 
<0  <0 
xs  p 

P 

to 

©  <0 
p  J3 
(0 

O  T> 

■H  C 
•O  <0 

c 

tI  h  • 

©  to 

MAH 

UBH 

C  <0  <0 
O  £2  £ 
fl  O  , 
a?  P 
H4J  O 
P  W  P 
to  3  o 
(0  U  P 
tv  xi  •a 
n  p  c 
o 

©  ©  o 
xi  £ 
o  p  © 
p  x: 
ft  c  P 
i  -h 
CD  XJ  XJ 
P  P 
g  -H  *0 
O  !t 
U  O 
<H  TJ  P 
Q) 

10  C  © 
P-H  O 

P  <0  M 
O  P  O 
P  ©  «W 
C  P 

o  in 

°Sx 

P  <0 
C  W-n 
0)  P 
P  <0  <W 
MHO 
P 

O  to  to 

•H  P 

•g  c 

<0  P  © 

to  c  c 
o  © 

rH  P 
O  P 
C  3 


O 

g 

o 


r.n  r\  n 


© 

P 

& 

•H 

IP 


6 


RDA-TR-14 42 00-002 


SEPTEMBER  1988 


A  computer-controlled,  time-  and  spatially-resolved  spectro- 
graphic  system,  shown  schematically  in  Figure  5,  was  used  for  ion 
temperature  and  electron  density  measurements.  In  order  to  make 
measurements  within  the  arcjet  thrust  chamber,  a  slit  1.6  mm  wide 
was  machined  in  the  outer  conductor  at  a  position  22.7  mm  from 
the  boron  nitride  backplate  that  allowed  observation  between  the 
cathode  and  anode.  To  minimize  the  perturbing  effects  on  the 
current  distribution  by  the  slit  a  series  of  4  identical  slits 
were  made  symmetrically  in  the  outer  conductor.  These  slits  are 
shown  schematically  in  Figure  6.  Time  integrated  line  profile 
measurements  were  made  of  the  argon  ion  434.803  nm  line  and  the 
696.28  nm  line  of  hydrogen  (HQ)  . 

If  the  velocity  component  of  a  radiating  particle  parallel 
to  the  direction  of  observation  is  vs,  then  the  wavelength 
appears  to  be  shifted  due  to  the  Doppler  effect  by  the  amount 

(1) 

where  aq  is  the  unshifted  wavelength  [2].  If  the  particle 
velocities  of  a  radiating  system  have  an  equilibrium  thermal 
distribution,  then  Doppler  broadening  results  in  a  Gaussian  line 
shape  and  the  observed  intensity  distribution  is  Gaussian  shape 
given  by 

1  <iX>  <2> 

where  aAq  =  ^aq;  v  =  (2kT/>)*  is  the  most  probable  velocity  of 
the  particles,  and  It  represents  the  total  line  intensity.  The 
full  width  at  the  half-maximum  intensity  point  (FWHM)  is  given  by 

AAj  —  2j m2  AAq  —  1.665  AAp 

AAj/Aq  =  7.713  X  10”5  (3) 

where  T  is  in  eV  and  n  is  the  reduced  mass  [2,3]. 

The  time  integrated  line  profile  of  434.803  nm  argon  ion 
observed  through  the  slit  in  the  outer  conductor  was  digitally 
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Figure  6.  Spectroscopic  measurements  were  made  internal  to  arc jet 
through  a  1.6  mm  slit  in  the  outer  conductor  at  the 
midplane.  Four  symmetrical  slits  were  cut  in  the  outer 
conductor  to  minimize  slit  effect  on  current  distribu¬ 
tion.  (This  figure  depicts  arcjet  with  cathode  and 
midplane  injection  of  mass  flow) . 
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recorded  with  the  spectrograph  in  second  order  vJth  the  entrance 
slit  at  30  microns.  The  system  resolution  was  0.002216  nm/pixel 
with  a  measured  instrument  profile  that  was  Gaussian  with  a  FWHM 
of  0.01262  ±  0.00045  nm.  The  digital  data  is  intensity 
integrated  over  the  line  of  sight  versus  wavelength  and  versus 
spatial  position  with  spatial  resolution  of  267  microns.  The  raw 
data  was  filtered  versus  wavelength  with  a  low  frequency  bandpass 
filter  with  frequency  cutoff  determined  by  the  allowable  spectral 
frequencies  from  the  instrument  profile.  The  data  was  then 
filtered  spatially  with  a  frequency  cutoff  of  0.35  cycles/mm. 

The  observed  line  profile  was  then  Abel-inverted  to  yield  local 
radial  line  profiles.  These  line  profiles  were  then  fitted  to  a 
Gaussian  and  the  FWHM  calculated.  Figure  7  shows  the  relative 
integrated  intensity  versus  radial  position  from  the  Abel 
inversion.  Figure  8  shows  the  Gaussian  fit  at  13.6  mm  from  the 
cathode  with  a  FWHM  of  0.02274  nm  with  an  error  of  +  0.00066  nm. 
Correction  for  instrument  profile  yields  a  FWHM  of  0.0189  nm 
which  corresponds  to  an  argon  ion  temperature  of  12.7  eV  +0.9  eV. 
Figure  9  shows  the  resultant  corrected  FWHM  versus  distance  from 
the  cathode  and  Figure  10  indicates  the  inferred  ion  temperature 
versus  distance  from  the  cathode.  These  inferred  ion  tempera¬ 
tures  are  very  high  and  are  approximately  10  times  the  electron 
temperature . 

Line  profile  measurements  were  also  made  of  hydrogen 
696.28  nm  line.  Hydrogen  comprised  1.5  percent  of  the  species 
injected  and  at  the  operating  conditions  interior  to  the  arc jet 
the  hydrogen  is  over  90  percent  ionized.  Thus,  the  observed 
hydrogen  line  intensity  was  very  low  and  required  the  sum  of 
30  shots  to  yield  reasonable  statistics.  The  raw  data  was 
filtered  and  Abel-inverted  in  the  same  manner  as  the  argon  data. 
Figure  11  shows  the  relative  integrated  intensity  versus  radial 
position  from  the  Abel  inversion.  The  original  objective  was  to 
use  the  measured  argon  ion  temperature  and  assume  that  the 
hydrogen  was  in  thermal  equilibrium,  then  deconvolve  the  hydrogen 
profile  using  the  calculated  Doppler  contribution.  The  decon¬ 
volved  profile  would  then  be  fitted  to  a  Lorentzian  profile  to 
yield  electron  density  from  Stark  broadening  calculations.  The 
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40.0 


RADIAL  POSITION  (MM) 

Figure  7.  Relative  intensity  of  argon  ion  emission  versus  radial 
position  from  Abel  inversion  of  Aril  434.803  nm 
emission  is  peaked  toward  the  cathode  (at  19.05  mm)  and 
monotonically  decreases. 


Figure  8.  Line  profile  of  Ar  II  434.803  nm  emission  at  13.6  ram 
from  the  cathode  shows  a  reasonable  fit  to  a  Gaussian 
profile . 
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Figure  9.  Calculated  FWHM  versus  distance  from  cathode  from 

argon  ion  emission  profiles  fitted  to  a  Gaussian  line 
profile  and  corrected  for  instrument  profile. 


profiles. 
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technique  was  not  possible  since  most  of  the  line  profile  was 
dominated  by  Gaussian  components  but  with  calculated  FWHM  much 
less  than  necessary  to  equal  the  argon  ion  "temperature".  The 
observed  Abel-inverted  hydrogen  profile  was  therefore  fitted  to 
a  Gaussian  and  to  a  Lorentzian  simultaneously.  The  accuracy  of 
this  method  is  not  as  good  as  if  the  two  broadening  mechanisms 
could  be  separated.  However,  regions  existed  where  there  was  no 
apparent  Lorentzian  contribution  and  had  a  good  fit  to  a  Gaussian 
thus  allowing  a  method  to  start  the  dual  fit.  Figure  12  shows 
the  calculated  hydrogen  temperature  from  a  fit  to  the  Gaussian 
contribution.  These  temperatures  are  certainly  more  reasonable 
than  calculated  from  the  argon  ion  line  profiles  with  a  peak 
calculated  temperature  of  only  0.85  eV.  A  fit  to  a  Lorentzian 
component  in  the  hydrogen  profile  only  occurred  within  10  mm  of 
the  anode.  The  resultant  calculated  electron  density  is  shown  in 
Figure  13.  Note  that  significant  errors  in  the  calculation  of 
the  Lorentzian  contribution  (i.e.,  50  percent)  can  be  tolerated 
and  still  yield  a  reasonable  measure  (factor  of  2)  of  the 
electron  density. 

If  the  particle  density  inside  the  channel  is  approximately 
2  X  1015/cm3  and  the  electron  temperature  is  on  the  order  of 
1.2  -  1.5  eV  then  the  hydrogen  and  argon  are  over  90  percent 
ionized.  The  coulomb  collision  frequency  of  hydrogen  ions 
(protons)  and  argon  ions  yields  relaxation  times  for  thermal 
equilibrium  of  tens  of  nanoseconds.  The  ion-neutral  relaxation 
times  for  the  same  species  are  significantly  longer,  on  the  order 
of  100  -  200  ns  but  are  still  less  than  the  expected  residence 
time  of  flow  particles  in  the  thrust  chamber.  One  therefore 
expects  better  coupling  between  the  hydrogen  atoms  and  the  argon 
ions  than  indicated  in  Figures  10  and  12.  The  calculation  of 
atom  temperature  assumes  a  Maxwellian  distribution  of  velocities 
characteristic  of  that  temperature.  Doppler  broadening  may  not 
be  necessarily  a  consequence  of  random  particle  motion,  but  can 
also  occur  as  a  consequence  of  collective  motion  in  separate 
regions  of  the  plasma  sampled  in  the  line  of  sight.  In  this 
case,  since  the  spectrograph  accepts  light  from  different  regions 
of  the  plasma  that  may  be  moving  in  complicated  ways,  the 
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Figure  11.  Relative  intensity  of  hydrogen  emission  versus  distance 
from  the  cathode  from  Abel  inversion  of  Ha  696.28  nm 
emission. 
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Figure  12.  Calculated  hydrogen  temperature  versus  distance  from 

the  cathode  from  Gaussian  component  of  line  profile  fit 
to  hydrogen  696.28  nm  emission. 
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Figure  13. 
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Figure  14 
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Calculated  electron  density  versus  distance  from  the 
cathode  from  Lorentzian  component  of  line  profile  fit 
to  hydrogen  696.28  nm  emission. 


DISTANCE  FROM  CATHOOE  (am) 

.  Average  relative  macroscopic  velocity  component  along 
the  line  of  sight  from  FWHM  measurements  of  Gaussian 
line  profiles  of  argon  ion  (+)  and  hydrogen  atom  (0) 
emissions  with  no  assumed  effects  due  to  kinetic 
temperature . 
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analysis  of  the  line  shape  cannot  show  whether  the  broadening  is 
a  consequence  of  the  thermal  motion  of  the  ions  or  atoms  or 
whether  it  is  due  to  a  distribution  of  radial  flow  (including 
fluid  turbulence) .  There  is  the  possibility  that  the  observed 
Doppler  broadening  of  the  argon  and  hydrogen  is  due  to  unresolved 
macroscopic  motion  which  may  result  in  Gaussian  profiles.  In 
this  case,  the  average  relative  macroscopic  velocity  component 
along  the  line  of  sight  is  related  the  line  profile  width  by 

c  aa  j 

<V*>av  =  TjTrt  T^~ 

and  has  nothing  to  do  with  the  kinetic  temperature  in  the  formula 
for  the  usual  Doppler  width  [4].  If  this  is  assumed  to  be  the 
case  for  the  observed  hydrogen  and  argon  profiles  then  the 
average  relative  velocities  would  be  equal.  Figure  14  is  a  plot 
of  velocity  versus  distance  from  the  cathode  for  hydrogen  and 
argon  using  equation  4)  and  assuming  no  contribution  to  the  line 
profile  due  to  effects  of  kinetic  temperature.  The  rough 
agreement  of  average  relative  velocities  between  hydrogen  and 
argon  tends  to  imply  that  the  mid-channel  kinetic  temperature 
is  low  with  the  presence  of  significant  macroscopic  motion  or 
fluid  turbulence. 

Examination  of  the  boron  nitride  backplate  after  many  shots 
indicated  that  cathode  material  was  being  deposited  on  the 
insulator  surface.  Such  depositions  could  contribute  to  erratic 
arc  behavior  and  might  also  distort  the  basic  current  distri¬ 
bution  within  the  thrust  chamber.  The  gas  inlet  through  the 
boron  nitride  backplate  was  therefore  redesigned,  as  shown  in 
Figure  15,  to  provide  more  uniform  gas  injection.  Using  the  same 
operating  conditions  as  with  the  previous  design  (i.e.,  8  kV 
charge  on  the  PFN  and  6  g/s)  the  line  profile  measurements  of  the 
argon  ion  at  434.803  nm  were  repeated  with  observation  through 
the  same  mid-channel  slit.  Figures  16a,  b,  and  c  are  contour 
plots  of  equal  intensity  versus  wavelength  versus  radial  position 
for  a)  the  raw  data,  b)  smoothed  data,  and  c)  Abel  inverted 
data.  The  smoothed  data  clearly  show  line  broadening  near  the 
cathode  but  the  Abel  inversion  lost  the  true  line  profile. 
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Figure  15.  Redesign  of  gas  feed  for  MPD  arc jet  with  gas  injection 
from  16  nozzles  at  mid-radius  through  backplate  and 
equal  gas  injection  from  a  channel  at  cathode  insulator 
interface  region. 


Figure  16.  Equal  intensity  contour  naps  of  Aril  434.803  nn  emission 

observed  through  slit  in  outer  conductor.  The  cathode  surface 
is  on  the  left,  the  anode  surface  on  the  right,  with  2.17 
mm/div  horizontal  and  0.004434  nm/div  vertical.  Figure  a)  raw 
data,  b)  smoothed  data  (note  that  peak  intensity  is  between 
10.8  -  15.2  mm  from  the  cathode),  c)  Abel  inverted  data. 
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Figure  17  shows  the  wavelength  integrated  relative  intensity 
versus  distance  from  the  cathode.  This  result  differs 
significantly  from  the  result  with  the  previous  design  shown  in 
Figure  7.  Again  fitting  the  observed  line  profile  to  a  Gaussian 
and  assuming  a  Maxwellian  velocity  distribution  yields  the  ion 
temperature  versus  distance  from  the  cathode  as  shown  in 
Figure  18.  The  dotted  line  in  Figure  18  at  0.3  eV  represents  the 
estimated  minimum  detectable  ion  temperature.  No  detectable  ion 
temperature  was  observed  between  7  -  19  mm  from  the  cathode  (note 
that  the  electrode  gap  is  23  mm) .  In  a  region  within  3  mm  of  the 
anode  (outer  conductor)  the  ion  temperature  rises  to  approximate¬ 
ly  2  eV.  Near  the  cathode,  the  ion  temperature  rises  to  above  2 
eV  by  6  mm,  but  due  to  problems  with  the  Abel  inversion  the 
accuracy  in  this  region  is  low. 
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Figure  17.  Relative  intensity  versus  distance  from  the  cathode  from 
Abel  inverted  Aril  434.803  nm  emission  with  new  gas  inlet 
geometry . 
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Figure  18.  Calculated  ion  temperature  versus  distance  from  cathode 
from  Gaussian  fit  to  Abel  inverted  line  profile  of  Aril 
434.803  nm  emission.  The  dashed  line  at  0.3  eV  represents 
minimum  detectable  temperature. 
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IV.  CONCLUDING  REMARKS 

The  experimental  program  to  date  has  provided  strong 

I 

indications  that  the  flow  within  the  MPD  arc jet  Qhamber  can  be 
dramatically  affected  by  the  choice  of  inlet  geometry.  In 
particular,  it  appears  that  introduction  of  mass  flow  near  the 
cathode  base  can  offset  the  development  of  a  strong  radially- 
inward  flow.  Such  flow,  as  suggested  by  the  linewidth 
measurements  in  argon  and  hydrogen,  would  have  important 
consequences  for  thrust  efficiency  since  shock  and  viscous  losses 
near  the  cathode  surface  would  be  expected.  The  present 
measurements  are  comparisons  at  the  same  total  current  but  with 
an  abrupt  variation  of  mass  flow  distributions  by  changing  the 
geometry  of  the  mass  injection  arrangement.  It  may  be  expected 
that  the  relative  strength  of  the  radial  (vs  axial)  flow  can 
change  significantly  as  the  total  current  and  mass  flow  values 
are  varied  with  fixed  injection  geometry.  Such  change  in  flow 
direction  can  substantially  affect  the  MPD  thruster  performance 
for  reasons  that  cannot  be  inferred  from  simple  scaling  relation¬ 
ships  based  on  electromagnetic  thrust  and  total  mass  flow. 
Examination  of  MPD  performance  by  variation  of  operating  current 
and  mass  flow  rate,  (i.e.  external  parameters)  with  a  fixed 
device  geometry  will  thus  be  inaccurate  for  prediction  of 
potential  limitations.  Each  set  of  external  parameter  choices 
may  require  a  particular  arrangement  of  mass  injection,  electrode 
geometry  and  length  in  order  to  achieve  optimum  performance. 
Internal  flow  measurements  and  MACH 2  calculations  in  the  next 
year  should  provide  considerable  insight  for  determining  the  best 
combinations  of  thruster  geometry  and  operating  conditions. 
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